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The light propagation of a probe field pulse in a four-level double-lambda type system driven 
by laser fields that form a closed interaction loop is studied. Due to the finite frequency width of 
the probe pulse, a time-independent analysis relying on the multiphoton resonance assumption is 
insufficient. Thus we apply a Floquet decomposition of the equations of motion to solve the time- 
dependent problem beyond the multiphoton resonance condition. We find that the various Floquet 
components can be interpreted in terms of different scattering processes, and that the medium 
response oscillating in phase with the probe field in general is not phase-dependent. The phase 
dependence arises from a scattering of the coupling fields into the probe field mode at a frequency 
which in general differs from the probe field frequency. We thus conclude that in particular for 
short pulses with a large frequency width, inducing a closed loop interaction contour may not be 
advantageous, since otherwise the phase-dependent medium response may lead to a distortion of the 
pulse shape. Finally, using our time-dependent analysis, we demonstrate that both the closed-loop 
and the non-closed loop configuration allow for sub- and superluminal light propagation with small 
absorption or even gain. Further, we identify one of the coupling field Rabi frequencies as a control 
parameter that allows to conveniently switch between sub- and superluminal light propagation. 

PACS numbers: 42.50.Gy, 42.65.Sf, 42.65.An, 32.80.Wr 



I. INTRODUCTION 

The propagation of a light pulse through a dispersive 
medium has been extensively investigated 0,0,0,0]. Re- 
cently the amplitude and the phase control of the group 
velocity in a transparent media have attracted much at- 
tention. It is well known that the propagation velocity 
of a light pulse can be slowed down ja, |fj , can become 
greater than the light propagation velocity in vacuum, 
or even negative in transparent media 0,13- Superlu- 
minal light propagation is a phenomenon in which the 
group velocity of an optical pulse in a dispersive medium 
is greater than of the light in vacuum. Despite the name 
superluminal light, it is generally believed that no infor- 
mation can be sent faster than light speed c in vacuum 
as explained by Chiao Thus, a group velocity faster 
than c as reported here does not violate Einstein principle 
of special relativity. The superluminal light propagation 
has been investigated for many potential uses, not only 
as a tool for studying a very peculiar state of matter, 
but also for developing quantum computers, high speed 
optical switches and communication systems flOj . 

There have been only few experimental and theoretical 
studies which realized both superluminal and subluminal 
light propagation in a single system. Talukder et al. have 
shown femtosecond laser pulse propagation from super- 
luminal to subluminal velocities in an absorbing dye by 
changing the dye concentration [TlT ] . Shimizu ct al. were 
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also able to control a light pulse speed, with only a few 
cold atoms in a high-fincssc microcavity by detuning the 
laser frequency from a cavity resonance frequency-locked 
to the atomic transition 12]. Most schemes for a speed 
control in atomic systems involve changing the frequency, 
amplitude or phase difference of the applied fields. It 
was shown that switching from subluminal to superlumi- 
nal pulse propagation can be achieved via the intensity 
of coupling field [H Q 0, and via the relative 
phase between two weak probe fields 0] ■ The intensity 
control of group velocity has been attempted by using 
a lower level-coupling field in the three-level atomic sys- 
tem . In another study, a scheme based on four- level 
electromagnctically induced transparency for switching 
from subluminal to superluminal light propagation via 
relative phase between driving fields was introduced |l8j] . 
Recently we have used an incoherent pump field to con- 
trol the li ght propagation from subluminal to superlu- 
minal [lii. l20j|. Further, we have studied the propaga- 
tion of light pulses in three-level closed laser-interaction 
loop atomic systems, where we exploited the fact that 
changing the relative phase between applied fields can 
modify the absorption and dispersion properties of the 
medium |2l|. The double-A setup is another scheme 
which allows for a closed laser-field interaction loop and 
provides a very rich spectrum of phenomena based on 
atomic coherence. One reason for this are the various in- 
terfering excitation channels in such a system [22]. The 
double A system has been investi gate d in the content 
of amplification without inversion |23(, phase sensitive 
laser cooling [24| , the propagation of the pairs of optical 
pulses [2^, optical phase conjugation j2(|, phase control 
of photoionization [23 , phase control of electromagneti- 
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FIG. 1: (Color online) The four- level double- A type schemes 
considered in the analysis, (a) is driven by four laser fields 
indicated by the red arrows with Rabi frequencies gij (i £ 
{3, 4},j 6 {1,2}) that form a closed interaction loop. This 
loop gives rise to a dependence on the relative phase of the 
different laser fields, and in general makes a time-independent 
steady state of the system dynamics impossible, (b) is the 
corresponding level scheme where one of the driving fields 
has been removed ((732 = 0). Here, the laser fields do not 
form a closed loop. The spontaneous decays with rates 7^ 
(2 6 {1, 2} ,j £ {3, 4}) are denoted by the wiggly green lines. 



cally induced transparency |22j and coherent population 
trapping j^. Recently Morigi et al (2^ have compared 
the phase-dependent properties of the o (diamond) four 
level system with those of the double A system. The 
propagation o f lig ht in a double A medium was studied 
earlier in Ref [23. The authors considered the case of 
complete resonance and the coherent population trap- 
ping (CPT) condition initially fulfilled, and found the 
possibility of amplification of one of the frequency pairs. 
Korsunsky et al. present a theory of a continuous wave 
light propagation in a medium of atoms with a double A 
configuration |30j |. They have shown that, when the so- 
called multiphoton resonance condition is fulfilled, both 
absorptive and dispersive properties of a such a medium 
depends on the relative phase of the driving fields. 

A peculiarity of the closed laser interaction loop sys- 
tems is the fact that a certain initial atomic state is con- 
nected to another atomic state via several combinations 
of laser field interactions. On the one hand, this gives rise 
to a dependence of the optical properties of the system 
on the relative phase of the driving fields. This allows 



for a great control of the optical properties and has been 
exploited in various circumstances, as discussed above. 
The drawback, however, is that such a system in general 
does not have a time-independent steady state. A time- 
independent steady state is reached only if a particular 
linear combination of the detunings of all incident laser 
fields is zero, that is, if the so-called multiphoton reso- 
nance condition is fulfilled. This assumption was made 
in the previous studies. But for a meaningful definition 
of the group velocity of a probe pulse in a given medium, 
one has to require that the medium itself including the 
coupling fields, which in this context merely aid in the 
preparation of the medium, has properties independent 
of the probe pulse characteristics. While the probe pulse 
which is finite in time necessarily consists of different 
frequency components and thus interacts with the atom 
with different detunings, the coupling fields have to be 
kept fixed in the calculation of the probe field propaga- 
tion, since the different frequency components interact 
with the medium simultaneously. In particular, it is not 
possible to adjust one of the coupling field detunings such 
that the multiphoton resonance condition is maintained 
throughout the pulse propagation. Thus, it is impossi- 
ble to assume the multiphoton resonance condition in the 
evaluation of the propagation of a finite probe pulse, and 
a time-dependent analysis is required. 

Therefore in this paper, we investigate the light prop- 
agation in a four-level double lambda system as shown 
in Fig. ^a) both with and without assuming the multi- 
photon resonance condition. We apply laser fields to all 
four dipole-allowed transitions to generate a closed laser- 
interaction loop. In the time-independent case where the 
multiphoton resonance condition is fulfilled, we find sub- 
and supcrluminal light propagation, controlled simply by 
the relative phase of the driving fields. For the time- 
dependent study beyond the multiphoton resonance con- 
dition, we apply a Floquet decomposition to the in gen- 
eral time dependent equations of motion. By compar- 
ing the Floquet decomposition to the time-independent 
treatment, we identify the respective frequency compo- 
nents found in this expansion with different interaction 
pathways in the closed loop system. We find that the 
medium response oscillating in phase with the probe field 
in general is not phase-dependent . The phase-dependent 
process contributing to the probe field susceptibility only 
occurs at a specific frequency, and is thus likely to distort 
the shape of a short probe pulse considerably if this spe- 
cific frequency is contained in the wide frequency spec- 
trum of the short pulse. Thus we conclude that for pa- 
rameters violating the multiphoton resonance condition 
as in the probe pulse propagation, inducing a closed loop 
interaction contour may not be advantageous. Neverthe- 
less, we demonstrate that also under conditions without a 
closed loop as shown in Fig.^b), both sub- and superlu- 
minal light propagation is possible with gain. Finally, we 
identify the Rabi frequency of one of the coupling fields 
as a convenient control parameter which allows to switch 
between sub- and supcrluminal light propagation. 
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II. THEORETICAL ANALYSIS 
A. The Model 

We consider a four-level atomic systems in a double- 
A configuration as depicted in Figure ^ The scheme 
consists of two metastable lower states |1), |2) and two 
excited states |3), |4). The electric-dipole allowed tran- 
sitions |1) - |3), |2) - |3) and |2) - |4) are driven by 
three coherent laser fields, while a weak tunable coherent 
probe field with central frequency lu p = W41 is applied 
to the dipole- allowed transition |1) — |4). The sponta- 
neous decay rates from level \i) (i G {3,4}) to the levels 
\j) (j G {1, 2}) are denoted by 27^. The electromagnetic 
driving fields applied to transition \i) — \j) can be written 
as 



E i:j = E ij e ij e- i ^ t - ki ^ + ^ + c.c. 



(1) 



with amplitude Eij, unit polarization vectors ey, fre- 
quencies u>ij, wave vectors fey and absolute phase 0^. 
Note that we do not assume absolute phase control in 
our calculations, but only relative phase control between 
the different driving fields. 

The Hamiltonian in dipole and rotating wave approx- 
imation is given by |3lL l32j| 



4 2 



(H9ime- iaim \l){m\+ h. 



(2) 



Z=3 m=l 



where the corresponding Rabi frequencies are denoted by 
gij = E^ (e^ dij )/H with dy as the atomic dipole moment 
of the corresponding transition. The energies of the in- 
volved states are denoted £7, (i G {1,...,4}), and we 
define the transition frequencies ujij = (Ei — Ej)/K. The 
arguments of the exponential functions are given by 



a 



<j 



LOijt — kijf + 



'ij 



(3) 



In a suitable reference frame, the Hamiltonian can be 
written as 



V =h(A 32 - A 31 )p 22 - hA 31 p 33 
+ h(A 32 - A31 - A 42 )p 4 4 
- ft (531/031 + .932/532 + 342P42 
^ \~ h.C.1 . 



(4) 



Here, we have defined pij = \i)(j\ and the corresponding 
operator in the new reference frame is denoted by pij 
(i,j G {1, ...,4}). It is interesting to note that now 
the residual time dependence along with the laser field 
phases appears only together with the probe field Rabi 



frequency gn in the parameter $ given by 

$ = At - Kr + O , (5a) 
A = (A 32 + A41) - (A31 + A 42 ) , (5b) 

K = (k 3 2 + ku) - (fcsi + £42) , (5c) 

00 = (032 + 041 ) - (031 + 042) ■ (5d) 

The parameters A, K and 0o are known as the multi- 
photon resonance detuning, wave vector mismatch and 
initial phase difference, respectively. Note that in gen- 
eral it is not possible to find a reference frame where 
the explicit time dependence due to A vanishes from the 
Hamiltonian, such that for A ^ Ono stationary long-time 
limit can be expected. Using the notations Eqs. (JjjJ and 
Eq. ©, the transformation of the operator pn, which 
will be of particular interest later on, to the chosen in- 
teraction picture can be written as 



p 41 = e-** e"*" p4i 



(6) 



We further define pn as the coherence pn in a reference 
frame oscillating in phase with the probe field En as 



p 41 = e Mil p 41 = e l *m 



(7) 

From the Hamiltonian Eq. Q , and including the spon- 
taneous decay in Born-Markov approximation, the den- 
sity matrix equations of motion can be derived as 

d . . . _ . . 



dt 



Pll = W 3 lP31 - *53lPl3 + WnPn e 



;'<!> 



ignPue '* + 27i 3 p 33 + 2714,944 



dt 



P22 = ig 32 P32 - «332P23 + «542P42 
- «542P24 + 2723/O33 + 2724P44 , 



g^P33 = ~miP31 ~ ig 32 P32 + «33lPl3 
+ «.932P23 - 273/O33 , 

d 

— P12 = l(A 32 - A 3 i)pi 2 + ig 31 P~32 - «ff32Pl3 

+ ifl , 4iP42e l * - igA2p~u - T^p^ , 

d 

-Q t Pi3 = -iA 31 p 13 + ig 31 (p33 - pn) - w 32 pi2 

+ j^iPise 1 * - T 13 p 13 , 

d 

— = i(A 32 - A 3 i - A 42 )/5i4 + ig\x^ (hi - Pn) 
at 

~ WI2P12 + ig 3 lP~34 - Ti4Pl4 1 

d 

q^P23 = ~iA 32 p 23 + ig 32 (p 33 - p 22 ) - ig 31 m 

+ «542P43 ~ r 23 P23 , 

d „ » • 

gj.P24 = -«A 42 p24 + «ff42(P44 - P22) - ifl&./^ie' 



«532P34 ~ r 24 P24 



8 . 



dt 



P3A = -i(A i2 - A 32 )p34 + «33lPl4 + ig32P2i 



*54lP3ie'* - igt 2 P32 ~ T34P34 , 



(8a) 

(8b) 

(8c) 

(8d) 

(8e) 
1) 

(8f) 

(8g) 
(8h) 

(8i) 
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pu = 1 - pu 



P22 — P33 



(8j) 



We have further defined jj = 71 j + -f2j , and our chosen 
level scheme implies 71 = 72 = 0. F^ — (2ji + 2'jj)/2 
(i 6 {1,2} and j £ {3,4}) are the damping rates of the 
coherences on transitions \i) — \ j). 

Due to the explicit time dependence of Eqs. © via 
the parameter A in $, it is clear that in general the 
system does not have a constant steady state solution. If, 
however, the so-called multiphoton resonance condition 
A = 0, K = is fulfilled, then the equations of motion (jHJ 
have constant coefficients and thus a stationary solution 
in the long-time limit can be found, which for a suitable 
choice of parameters depends on the constant relative 
phase 4>q. 



B. Linear susceptibility and group velocity 

The linear susceptibility of the weak probe field can be 
written as 0, 



2Nd. 



41 



Pa{u P ) , 



(9) 



where N is the atom number density in the medium and 
X = x' + *x"- The real and imaginary parts of x( w p) cor ~ 
respond to the dispersion and the absorption respectively. 
This expression for the susceptibility usually is applied to 
a single probe field frequency oj p of a narrow-bandwidth 
continuous-wave laser field. If the probe field consists of 
short pulses such that the individual probe pulses have 
a non-negligible frequency width, then Eq. © describes 
the propagation of an individual frequency component of 
the probe pulse. 

For a realistic example, we consider the sodium D\ 
transition, in which the transition rate, dipole moment 
and atom density are 7 = 2ir x 9.76 MHz, dm = 2.1 x 
10" 29 Cm and N = 1.3 x 10 12 cm -3 , respectively. The 
probe field with — O.OI7 is applied to the system, 
leading to (2 N da) / (e E 41 ) ~ 1. Eq. © refers to the 
part of the coherence p4i(w p ) oscillating at a frequency of 
the incident probe beam. This coincides with the trans- 
formed coherence pa introduced in Eq. (Q. Thus for our 
choice of example parameters, we obtain 



xK) ^ pa 



(10) 



Therefore throughout our numerical study, we will dis- 
cuss the transformed coherence /S41 in order to be study 
the light propagation in the medium. 

We also introduce the so-called group index, n g = c/v g 
where v a , the group velocity of the probe field, is given 

by Mm 



V 9 = 



1 + 27r X 'K) + 2ir£ ; x , {u Jp ) 



(11) 



This expression shows that, for a small imaginary part 
x"{oj p ) and positive steep dispersion, the group velocity is 



significantly reduced. On the other hand, strong negative 
dispersion can leads to an increase in the group velocity 
and even to a negative group velocity. 



C. Multi-photon resonance approximation 

In this section, we assume the multi-photon resonance 
condition A = 0, K = to be fulfilled. Then, the coef- 
ficients of the density matrix equations Eqs. ijSJ) do not 
have an explicit time dependence, and the system has a 
stationary steady state. In the particular case 

A31 = A 32 = A41 = A 42 = , (12a) 

7 = 7i4 = 724 = 713 = 723 , (12b) 

a simple analytical expression for the steady state solu- 
tion of p4i up to leading order in g±\ can be found: 



Pa 



-^7532 531542 



2D 



l53i| 2 + I532I 2 



1542 1 J |532 I 54ie~ 



2D 2 



*7l53i| 2 |53 2 | 2 |542| 2 5lie^ 
2L> 2 



(13) 



where 



£=l53i| 2 -|542| 2 + 7 2 (l53i| 2 + l532| 2 + |54 2 | 2 ) • (14) 

Transferring to the frame rotating in phase with the 
probe field using Eqs. © and (Q, we obtain 



pa 



-1-1932 531 542 e""' 

2D 

H 3 (l53l| 2 + |532| 2 + |542| 2 ) |532| 2 541 

2D 2 



*7l53i| 2 |532| 2 |542| 2 54ie 2 ^ 
2D 2 



(15) 



First, we note that the whole expression for p^i oscil- 
lates with same frequency as the probe field, as evidenced 
by the lack of an overall time-dependent factor in /)4i, 
and therefore contributes to the probe field susceptibil- 
ity. The three contributions to Eq. (|15fl admit a sim- 
ple interpretation in terms of the involved physical pro- 
cesses. The first component is proportional to 532 331 (742, 
and thus corresponds to a closed interaction loop involv- 
ing a sequence of all four dipole- allowed transitions, or a 
scattering of the driving field modes into the probe field 
mode. This round-trip depends on the relative multipho- 
ton phase 0o, which therefore occurs in this contribution. 
The second term in Eq. (|15fl proportional to 514 repre- 
sents a direct scattering of the probe field off of the probe 
transition, involving an excitation and a de-excitation of 
the probe transition. Since no closed interaction-contour 
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FIG. 2: (Color online) Interpretation of the different contributions to the probe field susceptibility in terms of transition 
pathways, (a) represents the interaction loop leading to a scattering of the driving fields into the probe field mode, (b) is the 
direct scattering of the probe field off of the probe transition, (c) shows a counter-rotating term. The solid red arrows indicate 
coupling field transitions, the dashed blue line is a probe field interaction. 



via different laser fields is involved, this contribution does 
not depend on the relative phase <f>Q. The last contribu- 
tion in Eq. (|15|) proportional to <? l4 depends on twice the 
relative field phase cf>o can be interpreted as a counter- 
rotating contribution. This preliminary interpretation of 
the individual contributions will become more transpar- 
ent in the time-dependent analysis without the assump- 
tion of multiphoton resonance in the following Sec. Ill Dl 



D. Beyond the multi-photon resonance 
approximation 

We now drop the multi-photon resonance condition 
and evaluate the equations of motion Eqs. ((SJ) for the 
general case of time-dependent coefficients. The Eqs. (JHJ 
can be written in a compact form 



^-R + E = MR, 
at 



(16) 



where 

R = (Plli Pl2, Pl3i PlA, p2\ip22-,p23i P24 , 

P31 , P~32 , P33 , P34 , PAl , PA2 , PA3 ) T (17) 

is a vector containing the density matrix elements, and 

E =(-27i4, 0, 0, -i&e***, 0, -2 724 , 

0, -igt 2 , 0, 0, 0, 0, ^4ie- ?A *, ig i2 , 0) T (18) 

is a vector independent of the density matrix elements 
which arises from eliminating one of the state popula- 
tions from the equations of motion via the trace condition 
Tr(p) = 1. Note that we have introduced the notation 



.941 = 541 e 



(19) 



The matrix M follows from Eqs. JSJ). Both the matrix 
M and the E can be separated into terms with different 
time dependence |34j . 



E = E + 34iEie 



-iAt 



#41 S -l e 



iAI 



M = M + g4iMie^ At + faM-ie*** . 



(20a) 
(20b) 



where Eo, E±i, Mq and M±\ are time independent. Using 
these definitions in Eq. (|16jl . we obtain 



— R + E + ff4ie- lAt E+i + g^e^E-i 



d - 
dt J 

(Mo + g 4 ie- lAt Mi + gl^M^R . (21) 



It is important to note that the time dependence of 
Eqs. JHJ only arises from the parameter A, such that the 
coefficients are periodic in time. According to Floquet's 
theorem [3j|, the solution R therefore has only contri- 
butions oscillating at harmonics of the detuning A. The 
higher-order harmonics in this frequency expansion are 
suppressed by powers of the probe field Rabi frequency 
f/4i relative to the other frequencies involved in the sys- 
tem. Since we are interested in the case where the probe 
field is weak, we truncate the Floquet expansion after the 
leading first order and obtain as ansatz for the solution 
R: 



R = R + gae-^Ri + g^e^R-x 



(22) 



By using Eq. in Eq. 121|) and equating the coefficients 
oscillating at different harmonics of A, we obtain the 
solutions for Rq and R±i as 

R = Mo-%, (23a) 
Ri = (A/ + iA)- 1 (E 1 -A/ 1 7?o), (23b) 
R_i = (M -iA)- 1 (T,- 1 -M-- L Ro). (23c) 

The medium response is determined by the 13th com- 
ponent p4i of R. Transferring back p^i to the operator 
pii in the frame oscillating in phase with the probe field 
using Eqs. JBJ and Q, we obtain 



Pai — [i?o]l3 



„i<5> 



54i[-Ri]i3 + gli[R-i]i3 e 



2/* 



(24) 



Here, [x]i3 denotes the 13th component of x. It can be 
seen that the contribution proportional to [i?i]i3 oscil- 
lates in phase with the probe beam, and thus contributes 
to the probe beam susceptibility independent of the fre- 
quency of the incident driving field. For A / 0, the two 
other contributions proportional to [i?o]i3 and [^-i]i3 
oscillate at different frequencies due to the time depen- 
dence of $ and thus do not contribute to the probe beam 
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susceptibility. However, even in the case A — > 0, the term 
[i?i]i3 is distinct in that the two other contributions in 
general propagate differently due to the wave vector mis- 
match K. 

Using the result Eq. (|24|l . we are now also in the po- 
sition to in detail understand the expression Eq. (|15fl we 
previously derived assuming A = and K = 0. For this, 
we explicitly evaluate [i?o]i3 5 [-Ri]i3 and [i?_i]i3 using 
the parameters Eq. (|12f> as well as A = and obtain 



[fi ] 13 e^= ~ ng32 £ 31g42 e-^" 



2D 



(25a) 



941 [R 



lj 13 



*7 3 (Iffsil 2 + I932I 2 + I942I 2 ) |332| 2 941 



2D 2 



(25b) 



941 [-R-i]i3 e 



2i$ _ — (|g31 1 2 |g32 | 2 |g42| 2 )g|l ^_2J?r+2^o 



2D 2 



(25c) 



Note that D is defined in Eq. Thus the time- 

dependent Floquet decomposition in the limit A — > 0, 
K — > yields exactly the result Eq. i|15|) of the time- 
independent analysis, as expected. The first part of 
Eq. 1)15(1 representing the scattering of the driving fields 
into the probe field mode arises from [Rq] 13 , as shown in 
Fig-EI a )- This contribution in general does not oscillate 
at the probe field frequency, but rather at the combi- 
nation frequency ^31+ W42 — W32 of the three driving 
fields. This frequency coincides with the probe field fre- 
quency only under multiphoton resonance. The contribu- 
tion proportional to [i?i]i3 shown in Fig.JIJb) is in phase 
with the probe field for all values of A, and is indepen- 
dent of the relative field phase. It represents the direct 
scattering of the probe field off of the probe field tran- 
sition. The third contribution proportional to [i?_i]i3 
can be interpreted as a counter-rotating term which in 
the Floquet expansion differs by 2A from the probe field 
frequency, and is depicted in Fig. Gl^c) . 

As an important result, we thus conclude that the 
phasc-dcpcndcncc of the loop-configuration studied here 
is restricted to the multiphoton resonance condition A = 
0, because it arises from the scattering of the coupling 
fields into the probe field mode. Furthermore, it can be 
seen that all contributions but the direct scattering ac- 
quire an additional dependence on the wave vector mis- 
match K, which influences the spatial emission pattern of 
these contributions. In general, only the direct scattering 
contribution can be detected in propagation direction of 
the probe beam regardless of the separation of detector 
and the scattering atoms. Finally, we note that the first 
expression [i?o]i3 is independent of A, whereas [i?i]i3 and 
[i?_i]i3 have complicated dependencies on A. 
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FIG. 3: (Color online) Real (blue solid line) and imagi- 
nary (red dashed line) parts of the atomic density matrix 
element pa in a reference frame oscillating in phase with 
the probe field as a function of the Raman detuning S12 = 
A31 — A32 = A41 — A42. The real part corresponds to the 
probe field dispersion, whereas the imaginary part describes 
the absorptive properties. The relative initial phase differ- 
ence between the four applied laser fields is cj>o — in (a), 
cj>a = 7r in (b) and cj>o = n/2 in (c). The common pa- 
rameters are 2713 = 2714 = 2723 = 2724 = 70, and thus 
Tia = r i4 = r 23 = r 24 = 70, T 12 = O.O, and T 34 = 2 7o . The 
detunings are A32 = A42 = 0, and the Rabi frequencies are 
chosen as 331 = (732 = 542 = O.670 and £41 = O.OI70. 



III. RESULTS AND DISCUSSION 

We now turn to the numerical study of our results from 
the master equations (jHJ. In a first step we assume that 
the multiphoton resonance condition A = 0, K = is 
fulfilled. Then the steady state solution for Eqs. (JSJ ex- 
ist and depends on the initial constant relative phase 4>o 
of the driving fields [see, e.g., Eq. (|15|l for the special 
case of parameters Eq. (1121) ] . In Fig. [21 we plot the ab- 
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-4-2 2 4 
£12/70 

FIG. 4: (Color online) Real (blue solid line) and imaginary 
(red dashed line) parts of the contribution to the Floquet 
decomposition representing the direct scattering of the probe 
field off of the probe transition, g^i [Ri]is [see Eq. (1241 ] . This 
figure depicts one of the processes contributing to the results 
in Fig. [3] and the parameters are the same as in this figure. 
Note that this particular contribution is independent of the 
phase (po, such that only one subfigure is shown, consistent 
with the interpretation in Sec. Ill Dl 
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sorption and dispersion spectrum versus the Raman de- 
tuning <5i2 = A31 — A32 = A41 — A42, for different ini- 
tial phases 4>q. The common parameters are chosen as 
2713 = 2714 = 2723 = 2724 = 70, where 70 is the to- 
tal decay rate of one of the upper transitions to a lower 
transition. Then, r i3 = r 14 = T 2 3 = T 2 4 = 70, Ti2 = 0, 
and T34 = 270. Further, A32 = A42 = 0, and the driv- 
ing fields have Rabi frequencies 513 = 523 = .924 = O.670, 
5i4 = O.OI70. The initial phase values in Fig.[3]are <fro = 
in (a), <fro = 7r in (b) and 4>q = ir/2 in (c). The slope of the 
dispersion (real part) is positive for <fro = 0, and accom- 
panied by gain (negative imaginary part). For (f>Q = n, 
however, the slope of the dispersion becomes negative to- 
gether with an absorption peak. For (f>Q = 7r/2, the curve 
shapes of real and imaginary part are interchanged. The 
results show the phase-dependence of the probe absorp- 
tion for a closed-loop configuration in good agreement 
with previous results as in Ref. |22|. 

We now evaluate the individual contributions to the 
Floquet decomposition of the density matrix, Eqs. (|24|l . 
using the same parameters as above in Fig. J2J). From 
Eq. H24f) . it is obvious that the numerically dominant con- 
tribution arises from [i?o]i3 cxp(i$), since the two other 
contributions are suppressed by one power of the weak 
driving field Rabi frequency gn. Consistent with this 
interpretation, the results for [i?o]i3 exp(i$) are almost 
identical to the ones shown in Fig. [31 with the difference 
given by the two other terms in Eq. (|24fl and the ne- 
glected higher-order terms of order 0(gf 1 ). Therefore, in 
the following, we only show the two difference contribu- 
tions <74i [i?i] 13 and g\ l [R-i] 13 cxp(2i<f>). Fig.0]shows the 
contribution <?4i[i?i]i3, which is the process contributing 
to the results in Fig.[3]depicted in Fig.|^b). As discussed 
in Sec. Ill Dl this process is insensitive to the phase </>o, 
such that only one subfigure is shown for all values of 4>q. 



FIG. 5: (Color online) Real (blue solid line) and imaginary 
(red dashed line) parts of the contribution to the Floquet 
decomposition representing a counter-rotating contribution, 
gli [_R_i]i3 exp(2i"I>) [see Eq. I|24|l ]. This figure depicts one 
of the processes contributing to the results in Fig. |S| and the 
parameters are the same as in this figure. This particular 
contribution depends on the phase 2<j)o, such subfigure (a) 
shows the result for cj>o — 0,7r and (b) for <j> — n/2, consistent 
with the interpretation in Sec. Ill Dl 



Note that the results confirm our expectation that the nu- 
merical contribution of this process to the full results in 
Fig. El is small. The contribution gl^R-i}^ exp(2i$) is 
shown in Fig. [SJ the corresponding process is depicted in 
Fig-IHc). This process depends on 2(j> , such that subfig- 
ure (a) shows the identical results for <j>o — and <fio = 71", 
whereas subfigure (b) shows the result for <fio = ir/2. 
Also the numerical contribution of this process to Fig. |3| 
is small as compared to [-Ro]i3 exp(i$). 

As discussed in Sees. [I] and III Bl the slope of the dis- 
persion spectrum against the probe field detuning has a 
major role in the calculation of the group velocity. But 
to define the group velocity of a probe pulse in a given 
medium properly, one has to require that the medium 
itself including the coupling fields has properties inde- 
pendent of the probe pulse characteristics. Note that in 
this context the coupling fields prepare the medium, and 
in that sense belong to the medium in which the probe 
pulse propagates, and thus their detunings needs to be 
kept fixed. The probe pulse, however, has a finite dura- 
tion and thus consists of different frequency components 
which interact with the atom with different detunings si- 
multaneously. In particular, it is not possible to, e.g., 
adjust one of the coupling field detunings such that the 
multiphoton resonance condition is maintained through- 
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FIG. 6: (Color online) Real (blue solid line) and imaginary 
(red dashed line) parts of the contribution to the Floquet de- 
composition representing the direct scattering of the probe 
field off of the probe transition, g±\ [J?i]i3 [see Eq. I24H 1. 
The susceptibility is plotted against the probe field detun- 
ing A41. This corresponds to a calculation relevant for the 
evaluation of the group velocity, but violates the multipho- 
ton resonance condition and thus requires a Floquet anal- 
ysis. The parameters in (a) are A31 = A32 = A42 = 0, 
2713 = 2714 = 2723 = 2724 = 70, 331 = I.870, 332 = O.270, 
342 = O.570, and 341 = O.OI70. Subfigure (b) shows the case 
A31 = IO70, A32 = A42 = 0, 2713 = 2714 = 2723 = 2724 = 70, 
fl3i = 332 = O.I70, 342 = O.570, and 341 = O.OI70. 
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FIG. 7: (Color online) Real (a) and imaginary (b) parts of 
the scaled contribution to the Floquet decomposition repre- 
senting the direct scattering of the probe field off of the probe 
transition, [.Ri]i3 [see Eq. 1241 1. The susceptibility is plotted 
against the probe field detuning A41 with other detunings 
chosen as A31 = A42 = 0. One of the Rabi frequencies is 
chosen as zero, gzi = 0, such that the setup corresponds 
to the one shown in Fig 0^b). The solid blue line is for 
331 = O.770, the red dashed line for (731 = O.8570, and the 
green dash-dotted line for 5131 = I.570. The other parame- 
ters are 2713 = 2714 = 2723 = 2724 = To, 342 = O.270, and 
341 = O.OI70. 



out the pulse propagation. Therefore, the propagation of 
a pulse cannot be described assuming the multiphoton 
resonance condition, and our numerical results Fig. 13151 
presented so far are not sufficient to evaluate the group 
velocity. 

Therefore, in Fig. [5^a). we show the probe field sus- 
ceptibility against the probe field detuning A41, while 
the coupling field detunings are A31 = A32 = A42 = 0. 
This corresponds to a calculation relevant for the eval- 
uation of the group velocity, but violates the multipho- 
ton resonance condition A = for most values of A41 
shown in the figure and thus requires the use of the time- 
dependent Floquet analysis. Consequently, in Fig.|^l we 
only show the component [i?i]i3 since this part of the 
Floquet decomposition is the only one which oscillates 
in phase with the probe field under these conditions. It 
can be seen in Fig. HJa) that around A41 = ±270, the 
real part of the susceptibility has positive slope, while the 
imaginary part is strongly negative. This indicates sub- 
luminal light propagation with gain. At about A41 = 0, 
the real part of the susceptibility has negative slope over 
a wide frequency range, together with a small positive 



or even negative imaginary part. This corresponds to 
superluminal propagation with small absorption or gain. 

If, however, the parameters and the laser field align- 
ments are such that for a particular value of A41 within 
the resonance structure in the susceptibility the multi- 
photon resonance condition is fulfilled, then the coupling 
fields may be scattered into the probe field mode. This 
would give rise to a huge, very narrow contribution to 
the probe field susceptibility [see Fig. |3] , which is likely 
to completely distort the probe pulse shape. In Fig.EJa), 
this could occur at A41 = if the laser field wavevectors 
are chosen suitably. Therefore for probe pulse propa- 
gation, parameters may be favorable which suppress the 
scattering of the coupling fields into the probe field mode. 
But under such conditions, the system is not phase- 
dependent. 

In Fig.[|Jb), we show a case where the multiphoton res- 
onance condition A = is violated for all shown values 
of A41 . This avoids the scattering of the coupling fields 
into the probe field mode, and is accomplished by choos- 
ing the detunings as A31 = 10 70, and A32 = A42 = 0. 
Around A41 = 0, subluminal light propagation with lit- 
tle absorption can be achieved. The fact that this and 
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similar result can be obtained for a rather large detuning 
A31 = 100<?3i as compared to the corresponding Rabi 
frequency then suggest that a closed-loop configuration 
of the laser fields may not be necessary at all. 

Therefore, in the final part, we discuss the dispersion 
and absorption spectrum for the system without a closed 
interaction loop shown in Fig. ^b) . Fig. [7] shows again 
the probe field susceptibility due to the direct scatter- 
ing of the probe beam off of the probe transitions, with 
Rabi frequency 932 = 0. Subfigurc (a) shows the real 
part, (b) the imaginary part. The solid blue line is for 
cjsi = O.770, the red dashed line for p 31 = 0.857o, and 
the green dash-dotted for g^i = I.570. All coupling 
field detunings are zero, and the other parameters are 
2713 = 2714 = 2723 = 2724 = 70, 342 = O.270, and 
,941 = O.OI70. If the Rabi frequency 531 for transition 
1) — |3) is small, then the results are similar to the sim- 
ple three-level electromagnctically induced transparency 
(EIT) case. The slope of the dispersion is positive, and 
subluminal light propagation appears with an EIT dip in 
the absorption spectrum, see the solid blue line in Fig.0 
With increasing Rabi frequency, the slope of the disper- 
sion around zero detuning becomes negative and super- 
luminal light propagation sets in. Together with super- 
luminal light propagation, zero absorption or even gain 
is achieved. Thus the intensity of the pump fields can 
be used as a simple control parameter to switch the light 
propagation from subluminal to superluminal. 

Fig. [8] shows another set of results for the system 
in Fig. Gib), with parameters A42 = — 570, A31 = 0, 
.932 = 0, 2713 = 2714 = 2723 = 2724 = 70, 342 = O.670, 
and f/41 = O.OI70, and 531 = O.770 for the solid blue line, 
331 = O.8570 for the red dashed line, and 931 = O.8570 for 
the green dash-dotted line. Due to the different coupling 
field detuning A42, the interesting probe susceptibility 
structure is shifted to A41 w — 570. In this region, for 
<7i3 = O.87, one finds subluminal light propagation ac- 
companied by strong gain. Increasing the Rabi frequency 
<7i3, again the light propagation switches to superluminal, 
still together with gain. 



IV. DISCUSSION AND SUMMARY 

We have discussed the propagation of a probe pulse 
through a medium that is driven by coupling and probe 
laser fields that form a closed interaction loop. This gives 
rise to a dependence of the system on the relative phase 
of the various laser fields, but in general prohibits the 
existence of a time-independent steady state of the sys- 
tem. A stationary steady state only exists if the so-called 
multiphoton resonance condition A = [see Eq. (JSJ] is 
fulfilled. A probe pulse, however, has a finite frequency 
width. Therefore, a time-independent analysis is insuffi- 
cient to correctly describe the pulse propagation through 
the medium. To solve this problem, we have solved the 
time-dependent system without assuming the multipho- 
ton resonance condition by means of a Floquet decom- 
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FIG. 8: (Color online) As in Fig. |7| but with parameters 
<?3i = O.770 for the solid blue line, 531 = O.8570 for the red 
dashed line, and g si = O.8570 for the green dash-dotted line. 
The other parameters are A42 = — 570, A31 = 0, (732 = 0, 
2713 = 2714 = 2723 = 2724 = 70, 342 = O.670, and #41 = 
O.OI70. The solid vertical line indicates A41 = 470. 



position of the equations of motion. We found that the 
different Floquet components can be interpreted in terms 
of different scattering processes. The phase dependence 
arises from a scattering of the coupling fields into the 
probe field mode, and thus only occurs at a specific probe 
field frequency. If this frequency is within the frequency 
width of the probe pulse, a strong distortion of the pulse 
shape can be expected. In some cases, it may be possible 
to find probe pulses that are sufficiently narrow in the 
frequency domain such that they fit into the frequency 
range where the scattering of the coupling fields into the 
probe field mode dominates. Especially for shorter and 
thus broader pulses, however, we conclude that it may be 
advantageous to avoid the scattering of the coupling fields 
into the probe field mode. Then, however, the system is 
no longer phase dependent, and a closed interaction loop 
may not be necessary at all. We expect that similar re- 
sults hold for other level schemes which involve a closed 
laser field interaction loop. 

Apart from these general considerations, using our 
time-dependent Floquet analysis, we have shown that 
for realistic parameter sets both the closed-loop double- 
A system and the corresponding system without a closed 
loop allow for sub- and superluminal light propagation 
with small absorption or even with gain. Further, we 
have identified the Rabi frequency of one of the coupling 
fields as a convenient control parameter to switch the 
light propagation between sub- and superluminal light 
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propagation. 
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